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Abstract. An example of exceptional points in the continuous spectrum of a 
real, pseudo-Hermitian Hamiltonian of von Neumann-Wigner type is presented and 
discussed. Remarkably, these exceptional points are associated with a double pole in 
the normalization factor of the Jost eigenfunctions normalized to unit flux at infinity. 
At the exceptional points, the two unnormalized Jost eigenfunctions are no longer 
linearly independent but coalesce to give rise to two Jordan cycles of generalized bound 
state eigenfunctions embedded in the continuum and a Jordan block representation of 
the Hamiltonian. The regular scattering eigenfunction vanishes at the exceptional 
point and the irregular scattering eigenfunction has a double pole at that point. In 
consequence, the time evolution of the regular scattering eigenfunction is unitary, 
while the time evolution of the irregular scattering eigenfunction is pseudounitary. 
The scattering matrix is a regular analytical function of the wave number k for all k 
including the exceptional points. 
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1. Introduction 

Non-Hermitian Hamiltonians are widely used to describe open quantum sistems in many 
fields of science HH7]. For continuum, PT— symmetric, non Hermitian Hamiltonians 
on an infinite line Bender et al. |8jj9] showed that the eigenvalue spectrum is purely 
real when the strength of the non-Hermiticity is small. Recently, Yogesh N et al. have 
investigated the signature of PT— symmetry breaking in coupled waveguides [10]. On 
the other hand, the readings in a measuring device are real numbers which, according 
to Quantum Mechanics, correspond to the points in the spectrum of an operator that 
represents an observable. This physical condition is stated in mathematical form by 
demanding that the spectrum of an operator representing an observable should be 
real m When the observable is represented by a self-adjoint operator the condition 
is automatically satisfied. But the reality of the spectrum of an operator does not 
necessarily mean that the operator is self-adjoint [121113] , This non-equivalence of 
selfadjointness and the reality of the spectrum of an operator was made evident by the 
discovery and subsequent discussion of a large class of non-Hermitian PT— symmetric 
Schrodinger Hamiltonian operators with a complex valued potential term but with real 
energy eigenvalues mmm- Even in the case of a radial Schrodinger Hamiltonian with 
a real potential term, the reality of the energy spectrum does not necessarily mean that 
the Hamiltonian is self-adjoint. In many cases, the spectrum of a real non-self-adjoint 
Hamiltonian is real but differs from the spectrum of self-adjoint ones in two essential 
features. These are, the possible presence of exceptional points and the possible presence 
of spectral singularities. 

Exceptional points already appear in the finite dimensional case of non-Hermitian 
Hamiltonian matrices depending on a set of control parameters [Tfil - i23] . whereas spectral 
singularities are characteristic features of Hamiltonians having a continuous energy 
spectrum [71I23H30]. Hence, they are not possible for finite dimensional operators. 
Exceptional points in the real continuous spectrum of a Schrodinger Hamiltonian with 
a real potential have received much less attention than in the finite dimensional case 
[6], [31], [32]. With the purpose of clarifying the topological nature of the exceptional 
points in the real, continuos energy spectrum of a quantum system, in the following 
we will present and discuss an example of exceptional points in a Hamiltonian with 
a real potential of von Neumann-Wigner type. This potential is generated from the 
eigenfunctions of a free particle Hamiltonian by means of a four times iterated Darboux 
transformation when the transformation functions are degenerated in the continuum of 
eigenfunctions of the free particle Hamiltonian [331135] . 

This paper is organized as follows: In section 2, we generated a Hamiltonian H [4] 
by means of a four times iterated and completely degenerated Darboux transformation 
which has two exceptional points in its real and continuous degenerated spectrum. In 
section 3, we compute the Jost solutions of H[ 4] normalized to unit probability flux at 
infinity. In section 4, we show that at k = ±q , the Wronskian of the two unnormalized 
Jost solutions of H[ 4] vanishes, this property identifies these points as exceptional points 
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in the spectrum of the Hamiltonian H[ 4], Section 5, is devoted to show that at the 
exceptional points, the Hamiltonian H[ 4] has a Jordan block matrix representation 
and a Jordan cycle of generalized eigenfunctions. In section 6 , we show that the regular 
scattering solution vanishes at exceptional points, while the irregular scattering solution 
has a double pole at k = ±q and the scattering matrix S(k) is a regular function of k. 
In section 7, we show that the presence of an exceptional point in the spectrum of H [4] 
does not alter the unitarity evolution of the regular time dependent wave function. The 
pseudounitary evolution of the irregular time dependent wave function is examined in 
section 8 . A summary of the main results and conclusions are presented in section 9. 


2. The potential V[4] 


A Hamiltonian that has two exceptional points in its real and continuous spectrum may 
be generated by means of a four times iterated and completely degenerated Darboux 
transformation [361137] , Thus, and according with Crum’s generalization of the Darboux 
theorem [38] . the function 

_ W 2 ( 0 , 0 , 0 , 0 2 0 , 0 3 0 , e ±ikr ) 

v ’ [1 " w^a^d^a^) • () 

is an eigenfunction of the radial Schrodinger equation with the potential 


H[4] = H 0 - 2^ In Hd(0, <9,0, 0, 2 0, 0, 3 0). ( 2 ) 

In these expressions, Wi(0,0,0, ...0 3 0) is the Wronskian of the transformation 
function 0 (g,r) and its first derivatives with respect to q, and W 2 ( 0 , 0 , 0 , ... 0 3 0 , e ±lkr ) 
is the Wronskian of the transformation function, 0(g, r), its first three derivatives with 
respect to q and e ±tkr which is an eigenfunction of the free particle radial Hamiltonian 
Hq with eigenvalue E = k 2 . The transformation function is 


0(9, r) = sin(gr + %)), (3) 

and 5,0 is shorthand for 00/0 q. The phase shift S(q) in the right side of eq.® is a 
smooth, odd function of the wave number q. 

The Wronskian Wi(0, 0,0, 0 2 0, 0 3 0) = Wi(q,r) is readily computed from ([3]) and 
the potential V[A] is obtained from Wi(q,r) and its derivatives with respect to r, 

V[4] = —2 1 (wpq^w^q.r) - W[ 2 (q,r)). (4) 

An explicit expression for W\{q,r) is the following 

Wi(g,r) = 16 (g 7) 4 - 12 (g 7) 2 + 8 (g 3 7 2 )(g 7 ) - 12 (g 2 7 i ) 2 
+ 24[(g 2 7i)(g7) + ( 97 ) 2 ] cos 2 6 + 3 sin 2 29 
+ [16(g7) 3 — 12<27 — 12 g 2 7 ! — 4 g 3 7 2 ] sin 26, (5) 


where 


6{r) = qr + S(q), 7 (r) = 0,6> = r + 70 , 
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7 o =d q 5(q), 7 i =dg6(q), 

72 = d^S(q). ( 6 ) 

For large values of r, the dominant term in the right hand side of eq.flSi) is (gy) 4 
which is positive and grows with r as r 4 . Hence, for large values of r, W\(q,r) is a 
positive and increasing function of r. However, if the phase 5(q) is left unconstrained, 
Wi(q,r) could take a negative value at the origin of the radial coordinate (r = 0), in 
which case it should vanish for some positive, non-vanishing value of r, giving rise to a 
singularity of the potential V[4] at that point. 

A necessary condition for the validity of the method of the Darboux transformation 
is that the potential generated should not have any singularities that are not present 
in the initial potential. In the case under consideration, this condition means that the 
Wronskian Wi(q,r) should not vanish for any positive value of r. Therefore, to avoid 
the appearence of singularities in V[4] at finite values of r, we will put the condition 

Wi(q, 0) > 0. (7) 


In explicit form, we have 




- 12 9 


,d 2 S(q)\ 


2l 


dq 2 J 

x cos 25(g) + 3 sin 2 25(g) 

+ H*rr 

x sin 25(g). 


12 „Ad) 


12 g 


dq 2 J 


4 g 


3 c(35(g)\ 


dq 3 J 


( 8 ) 


Despite its formidable appearence, this equation may readily be integrated. First, to 
simplify the notation, we define a new function t(q) as 

t(q) := tan 5(g), (9) 


then 


sin 25(g) = 


2*(g) 

1 + £ 2 (g) 


( 10 ) 


and 


( 11 ) 


form 


cos2<5( ‘ ,) = TTWy 

Written in terms of i(g), the condition of absence of singularities in V[4] takes the 


4(l + t 2 (g)) Wr(g,0) = (-t(g) + gi,(g)) 3(-t(g) + qt q (q)) + 6g% q (g) 

+ 2g \ qg {q) -3g 4 ^(-t(g)+ gf 9 (g))) , (12) 


in this expression, t q is shorthand for dt/dq. 

Now it is evident from (TT2T) that if t(q) satisfies 


- t(q) + qt,(q) = 13, 


(13) 
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the equation (fT21) becomes an identity and the condition (J7J) is satisfied provided that 

12/3 2 


W 1 (q,0) = 


(14) 


(1 +t\q)) 2 

Integrating we get 

t{q) = aq~P, (15) 

according to equation (jUJ) 

S(q ) = arctan(ag — (3), (16) 

in these expressions a and ft are free parameters. 

Once 5(q) is known as an explicit function of q, the functions 70 , 71 and 72 are 
obtained from its first, second and third derivatives 


7o = 


7i = 


a 


72 = 


dSjq) _ 

dq 1 + (aq — (3) 2 ’ 
d?8(q ) 2 a 2 (aq — j3) 

dq 2 
d 3 S(q) 


(1 + (aq-P) 2 r 

2a 3 (l — 3 (aq — (3 ) 2 ) 


dq 3 (1 + (aq — (3 ) 2 ) 3 

With the help of these expressions and (j5J) we get for Wi(q,r 

12 f3 2 24 paq 


(17) 


Wi (a r') — --- 1 - 

lW ’ J (1 + (aq - P) 2 ) 2 (1 + (aq - P) 2 ) 2 

12aq ((aq) 2 + P 2 - 1) . 

1 sm 2qr 


(cos 2 qr — 1 ) 


16 


(1 + (aq- P) 2 ) 2 
4 , 4 aq 


(qr) + 
3 (aq)' 


(1 + (aq- P) 2 ) 


(qr) 3 + 


(1 + (aq- P) 2 ) 


■ (qr) 


12 


(qr) + 


6 (aq) 2 

(1 + (aq- P) 2 ) 2 
2 , 2a< l 


24 


(qr) + 


2 , 2aq(l - P(aq - P)) 


(qr) 


16 (qr) 3 + 


(1 + (aq- P) 2 ) 2 
3 1 30 ^ / \ 2 


(qr) + 


(1 + (aq- P) 2 ) 
cos 2 (qr + 5(q)) 
3 (aq) 2 


(qr) 2 

(qr) 


(qr) 


(1 + (aq — P) 2 ) (1 + (aq-p) 2 ) 2 

1 — 6 (aq — P) 2 + (aq — /3 ) 4 


- 12 (qr) sin2 (qr + 5(q )) + 3 f| _ f 0)2)2 

2 A(aq - P)(l - (aq - P) 2 ) 1 

x sm 2 2 qr -\ ---t- - - sm 2 gr cos 2qr . (18) 

(1 + (aq-pyy 

From this expression, we verify that the Wronskian W\(q,r), for r = 0, behaves as 

12 P 2 


W\(q, 0) = 


(19) 


(1 + (aq- p) 2 ) 2 

The asymptotic behaviour of Wi(q,r) for large values of r is determined by the 
term (gy ) 4 which grows as r 4 and is the dominant term in the right hand side of eq. (l5|) , 
with this result and equation (J4|), we get 

sin 2 (gr + 5(q)) 


lim I/[4](r) fa 8 q- 


+ o(r 2 j 


( 20 ) 
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V[4](r) 



Figure 1. This graph shows the potential Vj4](r) for the values of the parameters 
a = 1, /? = 3 fixed at q = 1 


Then, V[4] is a potential of von Neumann-Wigner type [35] . The figure 1 shows the 
behaviour of V[4\ as function of r. 


3. The Jost eigenfunctions of H[ 4] 

The two linearly independent unnormalized Jost eigenfunctions of H[ 4], that belong to 
the energy eigenvalues E = k 2 and behave as outgoing and incoming waves for large 
values of r are given en eq.(l). Notice that all terms in the last column of the Wronskian 
W 2 (p,d q (/), ...d q p,e ±tkr ) are proportional to e ±lkr . 

Hence, the Jost solutions of eq. (JTJ) take the form 

r) = Wi(q, r) W± ^ ) r)e± ' fcr ’ (21) 

where the function w ± (k,r ) is the reduced Wronskian defined as 


w ± (k,r)e ±lkr = W 2 {(j),d q (j),...,d 3 q (l),e ±lkr ). 


( 22 ) 


From this definition, it follows that w ± (k,r) is a complex function of its arguments 

w ± (k, r) = u(k, r) ± iv(k 1 r). (23) 

A straightforward computation of the Wronskian 143(0, dqf, e ±lkr ) allowed us to 

find the following explicit expressions for the functions u(k,r ) and v(k,r ) 


,4 , a 2 12 , „4\ 2 aq) 0 i/i4 a 2 r 2 „4\ 


u(k, r ) = 12(A; 4 + Qq 2 k 2 + q 4 ) 
/3aq 


x 


(1 + ( ag - 0) 2 ) 2 

/ ___\ 3 


(1 + (aq- (3) 2 ) 2 
cos2gr + 16(k 2 — q 2 ) 2 (qr) A 
6 (aq) 2 


+ 24 (k A - 4q 2 k 2 - g 4 ) 


(1 + (aq- P) 2 ) 

3 (aq) 3 

(1 + (aq- P) 2 ) 2 


(qr) + 


(1 + (aq- P) 2 ) 2 


(qr) 


(qr) 


+ 24 


(, k 4 - 4 q 2 k 2 - g 4 ) 
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and 


x ( (gr) 2 + 


2 aq 


(qr) - 2 (k 4 - g 4 ) 


x 


2 (l + (ag-/3) 2 ) 

ac t) ( ac l P) cos 2(gr + 5(g)) + 16(A 4 — g 4 ) 


(1 + (aq- P) 2 ) 2 

3 , 3a< l 


-(qr) 2 + 


3 (ag)" 


* ( (,r) + (1 + {aq - PY) ' ' (1 +(aq-pyy 

— 12(A; 4 — 4g 2 & 2 — g 4 )(gr) sin2(gr + 5(g)) 

24(n- g ‘). (aqf 


(qr) 


x 


(1 + (aq — /3) 2 ) 2 
ag(l + (ag) 2 - /3 2 ) - 


12(A: 4 - 4q 2 k 2 - g 4 ) 


(1 + (aq- p) 2 ) 2 
+ 3 (k 4 + 6 q 2 k 2 + g 4 ) 


sin 2gr 


1 — 6(ag — (3) 2 + (aq — /3) 4 


(1 + (aq- P) 2 ) 2 
2(aq — P)(l - (aq - (3) 2 )- 


sin" 2 qr 


v(k, r) = 24 qk 


(1 + (aq- p) 2 ) 2 

P 

(1 + (aq - P) 2 ) 2 


sin 2 qr cos 2 qr 


(24) 


+ («g) (g + 5 At) 


+ 


(/3 2 -40ag-l)(A; 2 + g 2 ) 


24 qk(k 2 + q 2 ) 


x 


/3(/? 2 - 4/3ag + (ag) 2 + 1) 


+ 


(1 + (aq- (3) 2 ) 2 
6 Z 

3(ag) 2 


+ 96gfc 3 


^(ag) s 


(1 + (aq- P) 2 ) 2 J 

Q 

x cos2gr + 64qk(k 2 — q 2 ) (qr) 3 + - -—-—( qr ) 2 


1 + (ag - /3) 2 


(1 + (aq- P) 2 ) 2 
3 , 3ag 


(qr) — 24qk(k 2 + g )(gr) + 32 qk(k 2 — q ) 


x (gr) 3 + 


r (gr) 2 + 


3(ag)" 


1 + (ag - /3) 2 (1 + (ag - /3) 2 ) 2 

2Aqk(k 2 + g 2 )(gr) cos 2(gr + 5(g)) 

2 ag 


(gr) 


96g 3 & (gr) 2 + 


1 + (ag - /3) 


96 qk(k 2 — q 


2 24 v«<?) 2 ( a <?-/?) 


(1 + (aq- p) 2 ) 2 


M r ) 


(qr) sin 2(gr + 5(g)) 


+ 


12 qk(k 2 + q 2 ) 
48 qk(k 2 — g 2 


2 . 2 d ac l ~ Pf + ^P ac l ~ 1 


(1 + (aq- p) 2 ) 2 
(aq) 2 


(1 + (aq- P) 2 ) 2 


sin 2 gr 


1 — 6(ag — P) 2 + (ag — Py 


12 qk(k 2 + g 2 ) 

4(ag — P)(l — (aq- P) 2 ) 2 


(1 + (aq- P) 2 ) 2 
sin 2 2 qr . 


sin 2gr cos 2 qr 


(25) 


(1 + (aq-P) 2 ) 2 

For large values of r, the asymptotic behaviour of w ± (k,r) is dominated by the 
highest power of r. From eqs. (|24|) and (l25]l we get 
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w ± (k, r) m 16(k 2 — g 2 ) 2 (gr) 4 jl + 0(r 1 )|e ± * fcr (26) 

and from eq. ffTSl) we get 

W\(q,r) & lQ(qr) 4 [l + Ofr- 1 )], (27) 


hence, for large values of r 

/=*=(*,r) » (k 2 - q 2 ) 2 [ 1 + 0(r~ 1 )]e ±ikr . (28) 

The factor {k 2 — q 2 ) 2 is the flux of probability current at infinity of the unnormalized 
Jost solutions. 

Therefore, the Jost solutions of H [4] normalized to unit probability flux at infinity 


are 


F ± (k,r) 


fHKr) 

(k 2 — q 2 ) 2 


1 

( k 2 — q 2 ) 2 Wi(q, r ) 


k 2 * q 2 - (29) 


4. Exceptional points in the spectrum of H[ 4] 

Each pair of linearly independent Jost solutions belongs to a point = k 2 , with k 2 q 2 , 
in the spectrum of H [4]. At the point E q = q 2 , the two unnormalized Jost solutions 
coalesce to give rise to a Jordan chain of two generalized bound state eigenfunctions of 

m ra. 

The Wronskian of the unnormalized Jost solutions is readily computed from (12T|) 
and ([23]) 

W ( f + {k , r), f~{k , r)) = -2 ik(k + q) 4 (k - g) 4 , (30) 

At the points k = ±q, the Wronskian of the two unnormalized Jost solutions of H[ 4] 
vanishes, then the two unnormalized Jost solutions are no longer lineraly independent, 
coalesce in one bound state eigenfunction embeded in the continuum. This property 
identifies the points k = ±g as exceptional points in the spectrum of the Hamiltonian 

m- 


5. Jordan cycle of generalized eigenfunctions and Jordan block 
representation of H[ 4] 


The normalized Jost solutions of H [4], as functions of the wave number k , may be 
written as the sum of a singular and a regular part 


F ± (k, r) = 


+ 


-0s (g, r)e =R ' 5( ^ x|(<?, r)e =F * <5( ^ -0s ( — r)e ±lS( - g ^ 


+ 


(k - q) 2 (k - q) 

x±B( - q ’ r)e±m,> . + h p k , r ), 


+ 


{k + qf 


(k + q ) 

explicit expressions for the functions fjB(q,r) and y^(g, r) as functions of r, are 


Mq,r) = 


24 g 2 


Wi(q,r] 


(31) 


{—2g 2 y 2 cos 6 + [gy + g 2 7 x] sin 9 + siir 9 cos 9 } (32) 
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r) = Xb(q, r ) + Ho^ B (q, r), 

where 

O ~ 

Xb (q, r) = - 7 -[— 2 g 3 7 3 sin 8 — 3 g 2 7 2 cos 8 + 3 g 7 sin 3 8 

U ,(>/• r) 

— 72 g 3 sin 8 + 37 i 7 g 3 cos 0 + 3 sin 2 8 cos 0 ]. 

The functions 70,71 and 72 are given in equations (fI7|) 


l<A B (q,r)| 2 



( 33 ) 


(34) 


Figure 2. This graph shows the bound state eigenfunction tpB(q,r) as function of r 
computed for q = 1 and the values of the parameters a = 1 and (3 = 3 


l+s(q,r)l 2 



Figure 3. The generalized eigenfunction xb(q,v), as a function of r computed for 
q = 1 and the values of the parameters a = 1 and (3 = 3. 


In the figures 2 and 3 we show the graphical representation of the generalized 
eigenfunctions if B (qjr) and +s(g,r) as function or r. 
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The functions r ) an d Xb(<L r ) are the elements of a Jordan cycle of generalized 
eigenfunctions of the Hamiltonian that belong to the same point k — q in its spectrum, 
and satisfy the set of coupled equations 

H[A]f> B {q, r) = q 2 'f B (q, r) (35) 

and 


H [4]xb(<J, r) = q l XB(q, r) + 2 #b(?, r). 


(36) 


This Jordan cycle of generalized eigenfunctions is associated with a Jordan block 
representation of the Hamiltonian H[ 4], This property is made evident if the equations 
(|35li and (j36j) are written in matrix form 


tf[4]l 2 x2 


'I'st'J.i - ) = 'Hpiqt'ipMl- r). 


(37) 


where 

and 

*»<«> = (£ £)■ (39) 
From (137)1 . it is evident that the matrix "Hb(<?) is a matrix representation of the 
Hamiltonian H[ 4] in the two dimensional functional space spanned by the generalized 
eigenfunctions { ifsiq , r), XB(q, r)}. 

This space is a subset of the rigged Hilbert space of continuous, complex functions 
of the variables (q, r ) with continuous first and second derivatives with respect to r, 
with r in the semi-infinite straight line 0 < r < oo. Therefore, the two dimensional 
subspace of functions spanned by the generalized eigenfunctions {f>B{q, r), Xs(0V r )} is 
in the domain of H[ 4], 

The matrix T-Lsiq) a Jordan block of 2 x 2 [4IJ and can not be brought to diagonal 
form by means of a similarity transformation with a unitary matrix (423. 

The real non-symmetric matrix 1-LB^q) is r)— pseudo-Hermitian [T3j 

^b(s) = (40) 


with 

o)- (41) 

Hence, also the Hamiltonian operator H[ 4] when acting on the functional space 
spanned by the generalized eigenfunctions {ifsiq, r ), Xb(Qj r)} is also q— pseudo- 
Hermitian. 
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In this section it will be shown that the regular scattering solution tjj s (k,r) vanishes 
at k — ±q, and the irregular scattering solution ij; is (k,r) has a double pole at the 
exceptional points k = ±q. 

The regular scattering solution of H[ 4] is given by 


i(j s (k,r) = - F (k,r) - S(k)F + (k : r) , 

k^q, 

(42) 

where the scattering matrix S(k) is 



0(]] f~{k, 0) w~{k, 0) u{k,0)- 

1 J f+(k, 0) w+(k, 0) u(k, 0) + 

iv(k, 0) 
iv(k, 0) 

(43) 


From the Jost solutions eq. (12T1) . we get 


i/>s{k,r) 


i 1 11 

2 (F-g 2 ) 2 lFi(g,r) w+(k,0) 
w~(k 1 0)w + (k,r)e' lkr . 


w + (k, 0)w (k,r)e lkr 


(44) 


a graphical representation of the regular scattering solution is shown in figure 4. 


h s (r) 



Figure 4. The graph shows the regular scattering solution as function of r for the 
values of the parameters in V[4] and k = 1.5. 


In order to make explicit the properties of the scattering solutions as functions of 
k, for k close to the exceptional points at k = ±q, we will write the reduced Wronskian 
w ± (k,r) as an expansion in powers of (k — q ) 


4 

w ± (k, r ) = u(k, r ) ± iv(k, r) = e ^ wf(q, r)(k — q) e , (45) 

1=0 

explicit expressions for the coefficients wf(q,r), as functions of their arguments are 
given in the appendix A. 
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Then, the terms in the right hand side of eq. (1441) may also be expressed as an 
expansion in powers of (k — q) 

4 4 


ips(k, r) = 


2 ( k 2 — q 2 ) 2 Wi(q, r) w + (k, 0) 


£=0 m =0 


EE wj(q,0)Wm{q,r) 


x e 


—i(k—q)r _ - 


w t (q,0 )w+(q,r)e l{k q)r (k - q) 


— nf t+rrL l 


(46) 


In appendix B, we have shown that the first three terms in the summation in the right 
hand side of eq. ()46l) vanish. 

It follows that 

8 i 


i>s{k,r) = 


(■ k~q ) 


2 (k + q) 2 Wi(q, r ) w + (k, 0) 


EE 


re 


—i{k—q)r 


£=3 m =0 

w t-rn(l< 0)w+(q, r)e* (t -«) r l (k - q ) c ~ 3 + 0(k - q) 


(47) 


with the restriction wf(q, r) = 0, for £ > 5. 

Therefore, at the excepcional point, the regular scattering solution vanishes, 

i>s(q,r) = 0. (48) 


The energy eigenfunctions at the exceptional point are the generalized bound state 
eigenfunctions -0s(ov r ) and ys(q r ,r). 

Let us now turn our atention to the irregular scattering solution, 


As(k,r) 


F~(k, r) + S(k)F + (k, r) , 


k ± q- 


(49) 


(Mr) 



Figure 5. The graph shows the irregular scattering solution as function of r for the 
values of the parameters in V"[4] and k = 1.5. 


In this case, we have 


is(k,r ) 


11 11 
“ 2{k 2 -g 2 ) 2 ILi(g,r) w+(k,0) 
+ w~(k, 0)w + (k, r)e lkr . 


w + {k ) 0)w ( k,r)e lkr 


(50) 


Notice that, when the expression (1451) is substituted for w ± (k,r ) in (j50j) we obtain an 
expansion of the right hand side of (I5D]I in powers of (k — q) similar to (1461) . but now 
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the terms proportional to w^(q, 0)w±(q, r), w^(q,0)wf(q,r) and wf(q, 0)wf(q, r) are 
added. Hence, the singular terms in the expressions (l3Tf for the Jost solutions add 
instead of cancelling which makes t/y s (/e, r) singular at the exceptional points. Hence, 
the irregular scattering solution fj is (k,r), as a function of k, has a pole of second order 
at the exceptional points. 

Substitution of the expressions (145|) for w ± (k,r ) in (T5Q|) and multiplying the result 
times (k — q) 2 and taking the limit k —> q gives, 

lirn (k - q) 2 if is {k, r) = ip B {q, r). (51) 

k—tq 


Therefore, the bound state eigenfunction , if B (q,r) is the residue of second order of 
the irregular scattering solution fj is (k,r) at k — q. 

From equations (1241) and ()25l) for the funciones u{k , r) and v(k, r ) evaluated in r = 0 
we get 

u{k, 0) = ^ 1 -pyy [/3/c 4 + (6/3 - 20 aq)k 2 q 2 + (/3 - 4 aq)q A ] (52) 


and 


v(k, 0) = 


x 


48(3qk 

(1 + (ag-/3) 2 ) 2 

[(/3 2 — 4/3ag + 5 a 2 q 2 )k 2 + (/3 2 


A/3aq + a 2 q 2 )q 2 ] , 


(53) 


substitution of the equations (T52l) and (|53l) in (|43l) allows to write the following expression 
for the scattering matrix S(k ) 


S(k ) = exp2iA(k) (54) 

where the phase shift A {k) is 

4 qk [(/3 2 — 4(3aq + 5 a 2 q 2 )k 2 + (/3 2 — 4f3aq + a 2 q 2 )q 2 } 


A(k) = — arctan ■ 


/3/c 4 + (6/3 — 20 aq)q 2 k 2 + (/3 — 4 aq)q 4 


,(55) 



Figure 6. The phase shift A (fc) as function of the wave number k for q = 1 and the 
values of the parameters a = 1 and (3 = 3. 


the figure 6 shows the graphical representation of the phase shift. 

The scattering matrix is a regular analytical function of the wave number k, 
including the exceptional points, it does not have any singularity or pole. 
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7. Unitarity time evolution of the regular scattering eigenfunctions of H[ 4] 


The regular time dependent wave functions built as a linear combination of the regular 
scattering eigenfunctions of H[ 4] is 

T r (r, t) = j C{k)e~ l (r o) if s (k,r)dk. (56) 

This is a solution of the time dependent Schrodinger equation 

ih tttEldl = H{4}y r (r,t), (57) 

as it may be verified by substitution of (1561) in (1571) . 

The coefficient C\k) occurring in the right hand side of eq. (156D is a quadratically 
integrable function of k. 

Since the regular scattering eigenfunctions vanishes at the exceptional point 

*l>s{q,r) = 0, (58) 

it does not contribute to the integral in the right hand side of eq. d56lh 

Therefore, the presence of an exceptional point in the spectrum of H[ 4] does not 
alter the unitarity evolution of the regular time dependent wave function. 


8. Pseudounitarity time evolution of the irregular scattering eigenfunctions 


The two generalized eigenfunctions ijj B (q,r) and Xs(?; r ) belong to the same spectral 
point, E q = q 2 , in consequence, they evolve in time together. Hence, it should be 
convenient to introduce a matrix notation to deal with the two together, 


I'M) = r), 


(59) 


where Ts(g, r) is the two component vector of the doublet 


T B (g,r) 


f 4>B(q,r)\ 

\XB(q,r)J 


(60) 


and C(q, t ) is the 2x2 matrix of time dependent coefficent of the wave functions V ; b( , L r ) 
and XB(q,r). 

Substitution of T (r, t ) in the time dependent Schrodinger equation gives the 
following set of coupled equations written in matrix form 

ih^^-^ B (q, r) = C(q, t)H[4\l 2x2 V B (q, r) = C(q , t)H B (q)*B(q, r), (61) 
making abstraction of T B (g, r), we obtain 


where 


n B {q) 


q 2 0 
2 q q 2 


(63) 
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The matrix is a representation of the Hamiltonian H[ 4] in the 

two dimensional functional space spanned by the generalized eigenfunctions 
{Mq,r),XB(q,r)}. 

Integrating eq. (j62T) we get 

C(q,t) = (64) 

writing T-Lsiq ) in explicit form in (l64]h and computing the exponential, we obtain 

C(q,t) = e~ ^ • (65) 

Substitution of the expressions (1651) for C(q,t) and (l60li in (1591) gives 

ipB(r, t) = e~ l± ^f) B [q, r ) ( 66 ) 


and 


XB(r,t) = e 8 xs(q,r) 


• 2 q. 

i—te 

n 


iq t 




(67) 


9. Summary and conclusions 

An example of exceptional points in the continuous spectrum of a real, pseudo-Hermitian 
Hamiltonian H[ 4] of von Neumann-Wigner type is presented and discussed. The 
Hamiltonian H[ 4] and the free particle Hamiltonian H 0 are isospectral. In the general 
case, to each point in this continuous spectrum, correspond two linearly independent 
Jost solutions which behave at infinity as incoming and outgoing waves. Clearly in both 
cases, this continuous spectrum is doubly degenerate. However, here we have shown that 
in the continuous spectrum of H[ 4] there are two exceptional points at wave numbers 
k = ±q, these exceptional points are associated with a double pole in the normalization 
factor of the Jost eigenfunctions normalized to unit flux at infinity. At the exceptional 
points, the two unnormalized Jost eigenfunctions are no longer linearly independent 
and coalesce to give rise to Jordan cycles of lenght two of generalized quadratically 
integrablc bound states eigenfunctions embedded in the continuum and a Jordan block 
representation of the Hamiltonian H[ 4], The regular scattering eigenfunction vanishes at 
the exceptional point and the irregular scattering eigenfunction has a double pole with 
a second order residue equal to the bound state eigenfunction in the continuum and a 
Jordan bock representation of the Hamiltonian. In consequence, the time evolution of 
the regular scattering eigenfunction is unitary, while the time evolution of the irregular 
scattering eigenfunction is pseudounitary. The scattering matrix S(k) is a regular 
function of k at the exceptional points, that is, the Jordan cycle of generalized bound 
states eigenfunctions of H[ 4] in the continuum is not associated with a pole of the 
scattering matrix. 
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The dependence of w ± (k,r ) on k is readily made evident by expanding the reduced 
Wronskian defined in eq. (l45]i in powers of (k — q) 

4 

w ± (k,r ) = e Tie ^^ w f(q,r)(k — q) e , (A.l) 


e=o 


the coefficients of the powers of (k — q) in this expansion are 
Wq (q,r) = 4q 2 Wi{q,r)ip B (q,r) 

wf(q,r) =4qW 1 (q,r) ip B (q,r) + < 2 X B (< 2 ,r) = fiq^ B (q,r) 


(A. 2 ) 

(A.3) 


wf(q,r) = ~ 2Wi(q,r)f> B {q,r) cos 2 6 + 6qWi(q,r)xB(q,r) 

+ 16g 2 |4g 4 7 4 — 3g 2 y 2 — 3g 4 y 2 + 2g 4 yy 2 + 3 sin 2 9 cos 2 0 j cos 6 

T *21Ai(g,r)|(2g7 + g 2 7i)V’s(g,r) + 2g 2 7y s (g,r)| (A.4) 


uf(q,r) = 2W 1 (q,r)x B (q,r) - - Wi(q, r)if B {q, r) 

+ 8 q (( 8 g 4 7 4 — 18 g 2 7 2 + 3<2 3 7i7 — 6 g 4 7 2 ) cos 6 
+ ( 8 g 3 7 3 + 3 <27 - 12 g 2 7 i + 2 g 3 7 2 ) sin 3 9 

— 12 ( 2 g 2 7 2 + g 3 7 i 7 — 1)) sin 2 9 cos 9 

=F * 4 g 7 hhi(g, r)x B {q, r) + 2(7 + 7 i)H / i(< 2 , r)if B (q, r) 

— 8 q (( 12 g 3 7 3 — 6 g 2 7 i — 3 g 3 7 2 ) sin 2 0 cos 0 

+ 12 (g 2 7 2 + 3 g 3 7 i 7 ) cos 2 0 sin 0 — 12 g 3 7 3 cos 0 ) 


(A.5) 


cuy(g, r) = —-VFi(g, r)ip B (q, r) cos 0 + 8 ( 4 g 3 7 3 + g 3 7 2 ) sin 0 cos 0 


q 


— 48g 717 sin 0 


cos 0 ± i 


+ 16g 4 7 4 — 36g 2 7 2 + 12 g 2 7 i( 2 g 7 — g 2 7 i) + 36 sin 2 0 cos 2 0 

. 3 ^, „u ,2 £il _u „• —r)' 0 n(q r , r) cos 0 

- 0 2 

+ 8 ( 4 g 3 7 3 + g 3 7 2 ) sin 0 cos 0 + 16 g 4 7 4 — 36g 2 7 2 
+ 12<2 2 7!(2<27 — q 2 71 ) + 36 sin 2 0 cos 2 0 — 48<2 3 7i 7 sin 2 0 
x sin 0 


(A. 6 ) 


An expansion of the reduced Wronskian, w ± (k,r ), about the point k = — q is 
readily obtained from eq. (1A.1I) . From time reversal invariance of the Hamiltonian plus 
boundary conditions we get 

f + {~k,r) = f~(k,r) (A.7) 

and 

w^(k,r) — w T (—k,r). (A. 8 ) 
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Hence, from (lA.lj) 

4 

w ± (k, r ) = e Tl6 ^ l) l wf(q, r)(k + q) e , (A.9) 

e=o 

the coefficents wf(q,r) are explicitly given in eqs. (lA.2l - IA.6I) as entire functions of q. 
Their domain of analiticity is trivially extended to negative values of the argument q 
just by changing q by — q in (1A.2I - IA.6I) . Then it may be verihed that 

(-l) e wf(q,r) = wf(-q,r). (A.10) 

Therefore, 

4 

w ± (k, r ) = e ±te (~g) YX-iywt(\q\,r)(k + q) £ , k < 0. (A.11) 

£=0 


Appendix B. 


In this appendix it will be shown that the regular scattering eigenfunction vanishes at 
the exceptional points k = ±q. 

The regular scattering solution, eq. (l29j) . is 

1 1 


^s{k,r)= 1 -- 


w + (k,0)w (fc,r)e 


—ikr 


(B.l) 


2 (k 2 — q 2 ) 2 Wi(q,r) w+{k, 0) L 
— w~(k, 0)w + (k, r)e lkr . 

With the reduced Wronskian w ± {k ) r ) as an expansion in powers of ( k — q ), eq. flA. Ill , 
we get 

8 i 


Mk ’ r) = 5 W - 7? W,L r) „+(*. 0) 5 E K-<«'^ 
x 0) m +(9, j (k - q)>. 


(B.2) 


with the restriction wf(q,r) = 0 for £ > 5. 

Writing the first three terms in the summation, corresponding to £ = 0,1 and 2, 
and using the expansion of e -*(fc-9h, we 


fj s (k,r) = 


-{Co (q,r) 


2 (k 2 — q 2 ) 2 W\{q, r) w+(k, 0) 

+ Ci(g,r)(k -q) + C 2 (q,r)(k - q) 2 + 0{{k - q) 3 ) 

8 l 

+ J2 r)e~ i{k ~ q)r 


1=3 m =0 


(k~q)% 


- w e _ m (q,0)w+(q,r)e i{k q)r 
where the functions Co(? ; r ); Ci(?V r ) Y C 2 (<?, r) are given by 
Co(9, r) = - c.c. 


(B.3) 

(B.4) 


Ci( 9 .f) = -™J(?,0)tu o (j,r)r + ioJ'(j,0)io 0 (?,r) 
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+ wf(q,0)w 1 (q,r) - c.c. (B.5) 

C 2 (q, r) = - l;wf(q,0)wf(q, r)r 2 - iwf(q,0)wf(q, r)r 
— iwf (g, 0 )wf (g, r)r + wf (g, 0 )wf(q, r)r 

+ (q, 0)wf (g, r) + w£ (g, 0 )wf (g, r) - c.c. (B. 6 ) 

in these expressions, c.c. is shorthand for complex conjugate. 

From the eq. (1A.2I) , w^(q, r) = wf (q,r) = Wo(q,r), where wo(q,r) is the real 


function 

w 0 (q, r ) = 4g 2 tFi(g, r)ip B (q, r) (B.7) 

then, from eq. (1B.4D . 

Co(q,r) = 0. (B. 8 ) 

From eq. (1 A. 3h we note that the imaginary part of wf(q,r) is 

wf(q, r ) = =F4fFi(g, r)q 2 ^ B (q, r), (B.9) 

hence, (i(q,r) is written as 

Ci (q,r) = - 2iw 0 (q,0)w 0 (q,r)r - 8iq 2 'y o W 1 (q,0)'ip B (q,0)w o (q,r) 

+ 8 fg 2 7 w 0 (q, 0)Wi(g, r)ip B (q, r) (B.10) 

and using (IB. 71) we get 

Ci(g, r) = 2iwo(q, 0)tu 0 (g, r)(—r — 70 + 7 ). (B.ll) 

But, from eq.(| 6 j), 7 = r + 70 , then 

O(g,r) = 0. (B.12) 


The function ( 2 (g, r) is computed from the expressions for wf(q,r) and wf(q,r), 
we get 

C 2 (q,r) = - 8iqW 1 (q,0){'if B (q,0) + qxB(q,0)]w o (q,r)r 

- 8iqW 1 (q,r)[i/> B (q,r) + g%s(g, r)]w 0 (q, 0)r 

+ 32fg 2 fF 1 (g,O)['0i?(g,O) + gy B (g,O)]fFi(g,r)g 7 ' 0 B (g,r) 

- 32iq 2 W 1 (q,r)['f B (q,r) + qx B (q,r)}W 1 (q,0)qxo'f B (q,0) 

- 4:iWi(g, 0)[(2g7o + g 2 7i)tM<h 0) + 2g 2 7 0 x B (g, 0)]w 0 (g, r) 

+ 4 flFi(g,r)[( 2 g 7 + g 2 7 1 )V> i? (g,r)+ 2 g 2 7 x B (g,r)]u 7 0 (g, 0 ). (B.13) 
Using eq. (IB.71) for w 0 (q,r), we get 

C 2 (g,r) = 32fg 3 lUi(g, r)lUi(g,0) { 2 ^ B (g, r)-0 B (g, 0) 

+ g[V’B(g,r-)xB(g,o) + V’B(g,o)xB(g,r)]}(-r-7o + 7). (B.14) 

But, 7 = r + 70 , then 

C 2 (?,r) = 0. 


(B.15) 
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With these results, the regular scattering solution i/j s (k, r), as function of k, close to the 
exceptional point k — q, is 


A(k,r) 


1 1 k — q 1 

2 (k + q) 2 Wi(q, r) w + (k, 0) 


8 i 

EE wl_ m (q,0)w m (q,r) 

£=3 m =0 


x e-^ r - wf_Jq, 0) W +( g ,r) e ^> 
+ 0((A;-g)). 


(* - ?) £ - 3 


(B.16) 
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